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ACTIVATION ENERGIES FOR REACTIONS OF HYDROXYL RADICALS 
WITH HYDROGEN AND CARBON MONOXIDE 
by Edgar L. Wong and Frank E. Belles 
Lewis Research Center 


SUMMARY 


Direct measurements of the effect of temperature on the rates of reactions of hy- 
droxyl radicals are lacking in the literature. In the present work, rate constants for the 
k o 

H 2 + OH — =-*-H20 + H reaction were determined at temperatures of 304, 403, and 
503 K. A stirred reactor was used, and the reaction mixture was analyzed by means of 
a time -of -flight mass spectrometer. The results are described by the following 
Arrhenius expression: 
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k 2 = (4. 7±1. 2)xl0 ^ exp 


-(4200 ±300) 


RT 


cm^ mole * sec * 


The activation energy of 4. 2±0. 3 kilocalories per mole (17. 6±1. 3 kj/mole) was less 
than previous estimates, based on indirect evidence, of about 5. 15±0. 15 kilocalories per 
mole (21.56±0.63 kj/mole). . 

Earlier data obtained by the authors on the CO + OH =*- COg + H reaction were 

reanalyzed, using the new result to correct for the effects of the simultaneously occur- 
ring H 2 + OH reaction. The following, improved rate expression for the CO + OH re- 
action was obtained: 
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k t = (3. 5±1. 1)xl0 exp 


-(1120±80) 

RT 


3 , -1 -1 

cm mole sec 


INTRODUCTION 

The study of elementary reactions involving the hydroxyl radical (OH) is important 
for a better understanding of overall combustion processes. During the last few years, 




various investigators (refs. 1 to 5) have made quantitative OH reaction-rate measure- 
ments. With the exception of the work of reference 4 and a previous report by the 
authors (ref. 5), there have been no direct measurements of the effect of temperature on 
OH reactions. In reference 5, a mass-spectrometric stirred-reactor method was used 
to measure the rate of the following reaction at temperatures of 310, 440, and 610 K: 

k l 

CO + OH — U- C0 2 + H (I) 

Since the OH was generated by the reaction of partly dissociated Hg with nitrogen di- 
oxide (NOg), there was competition for OH because of the reaction 

k 2 

h 2 + oh — ^h 2 o + h (H) 
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This effect was accounted for by assuming a kg value of 2. 3x10 exp (-5200/RT) cubic 
centimeters per mole per second recommended in reference 2. However, the rate con- 
stants deduced for reaction (I) from our high -temperature data were not the same when 
two different methods of data reduction were used. This suggested that the activation 
energy recommended for reaction (H) might be in error and also made uncertain the acti- 
vation energy reported in reference 5 for reaction (I). 

The present work was therefore undertaken to measure directly the activation ener- 
gy for reaction (H) and to apply the result to a reanalysis of the data for the CO + OH 
reaction. 


EXPERIMENT AND PROCEDURE 
Apparatus 

The 300-cubic-centimeter stirred reactor and its connection to the time -of -flight 
mass spectrometer are shown in figure 1. A more detailed description may be found in 
reference 6. As in previous work, the inside surfaces of the reactor were coated with 
phosphoric acid. 


Gases 

Hydrogen, argon, and oxygen had stated purities of at least 99. 9 percent and were 
used without further purification. Nitrogen dioxide was about 90 percent pure and was 
further purified by fractionation at dry-ice temperature. 
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Calibration 




The rate of reaction (n) was determined from measurements of Hg and HgO concen- 
trations in the stirred reactor. Hydrogen (H 2 ) was metered, and it was only necessary 
to determine the sensitivity of the mass spectrometer to HgO. To guard against possible 
instrument drift, at least two calibrations were made before and after a run at a given 
temperature. 

Figure 2 is a typical calibration showing the checkpoints made before and after a 
kinetic rim. This water calibration was performed by saturating argon with distilled 
water at two temperatures: near 5° and near 15° C. The signal at zero water concen- 
tration is due to unavoidable background and was subtracted out of the results. In gen- 
eral, the mass-spectrometer water calibration curves obtained before and after a run at 
a given temperature agreed to within ±5 percent. 

The mass spectrometer was calibrated for Og by measuring the signals produced 
with metered well-mixed flows of this gas with argon. 


Production of Hydroxyl Radicals 

As in reference 5, OH was produced by titrating hydrogen atoms with nitrogen di- 
oxide (ref. 1) in the stirred reactor: 

H + N0 2 OH + NO (IH) 

Atomic hydrogen was made by subjecting a dilute (~2 percent) mixture of H 2 in argon to 
a microwave discharge. The amount of H 2 dissociated was observed to be 35 to 40 per- 
cent by monitoring the H 2 peak when the discharge was turned "off” and "on. " This in- 
formation was necessary for the following reasons: 

(1) In the present experiment it was assumed that the initial OH concentration |oh| 0 

is equal to the concentration of added NOg jNOgJ^ ' ^is assumption is necessary for cal- 
culating the steady-state concentration of OH [oiijt (to he discussed in the next section). 
To ensure that [ohJ 0 is equal to |n 0 2 J , the H atoms must always be in excess; if N0 2 

is in excess, there may be further hydroxyl formation resulting from titration of H atoms 

formed in subsequent processes, such as O + OH ► 0 2 + H. To be sure of an excess 

of H atoms, N0 2 was added in an amount equivalent to about one-half of the H atoms pro- 
duced in the discharge. 

(2) The information on the dissociation of H 2 was also needed to calculate the concen- 
tration of H 2 in the reactor when more Hg was added to study the Hg + OH ► HgO + H 

reaction. 
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The initial OH concentration was reduced during the time of residence in the reactor 
by the following two reactions, which, like reaction (HI), are also very fast (ref. 7): 


OH + OH — 

h 2 o + O 

(IV) 

O + OH — 

-o 2 + h 

(V) 


At room temperature and in the absence of any further additions of H 2 , reaction (II) is too 
slow to compete with reactions (IV) and (V) for OH. Under these conditions, the resi- 
dence time in the reactor was sufficient for reactions (in) to (V) to proceed close to com- 
pletion with the following overall stoichiometry (ref. 7): - 

2H + 3N0 2 3NO + H 2 0 + Og (VI) 

Some preliminary work on the H + N0 2 titration confirmed equation (VI); the result 
is shown in figure 3. It can be seen that the amounts of HgO and 0 2 measured by the mass 
spectrometer are accurately equal to one -third of the added NO 9 . (In fig. 3 as elsewhere 
in this report, flow rates reduced to standard conditions (F, cm^/sec, 300 K and 
750 torr) are used interchangeably with concentrations, since one can be changed to the 
other by means of the known total flow rate. ) 

Inasmuch as OH could not be measured directly because its concentration in the re- 
actor is too low, the result shown in figure 3 is important for two reasons: (1) it shows 
that the assumed equality of initial |ohJ o to [n 0 2 J is valid; and (2) it suggests that OH 

does not disappear, to any important extent, by ternary or wall reactions because such 
reactions would be unlikely to preserve the 1:3 ratios of both H 2 0:N0 2 and 0 2 :N0 2 called 
for by the gas-phase reaction (VI). 


Calculation of Rate Constants 

H 2 + OH reaction. - The rate of production of water in a stirred reactor in which re- 
actions (II) and (IV) are occurring is 

where At is the residence time in the reactor and the subscript t indicates the steady- 
state concentrations in the reactor. Equation (1) can be rearranged to obtain 
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( 2 ) 





V 


. /M 


At 



[ OH ]t[ H 2] t , 


In equation (2), the quantities ^Oj, and At are measured. However, as was 

stated earlier, joH ^ could not be measured and it was therefore necessary to calculate 
it in a manner similar to that used in reference 5. If it is assumed that atomic oxygen is 
in a steady state because of the great speed (ref. 7) of reaction (V), the rate of reaction 
of OH in a stirred reactor is given by 


4qhL (Ht-Ho) . 


At 


At 


3k 4 [o H ] t 2 + ^[OH],^ 


(3) 


By invoking equation (1), equation (3) can be rewritten as 



and 



(5) 


where [oh] q = [NOgJ^ This result was substituted into equation (2), which was then 
solved for kg . 

CO + OH reaction . - As was pointed out in reference 5, the rate constant for this 
reaction may be derived from the data of reference 5 in two ways. The first and more 
direct method follows directly from the reaction-rate expression for a stirred reactor: 



= k l[CO] t [OH] t 


( 6 ) 


The quantity [OH] t is calculated from an expression similar to equation (5) but with an ad- 
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ditional term accounting for the loss of OH in reaction (I): 


zNt At + 1 ) 2 + l2k 4 At ([ OH ]o - [ co : 
At ~~~ 




The second method makes use of the fact that each addition of CO was accompanied 
by a simultaneous decrease in the observed HgO peak intensity, because of the competi- 
tion for OH of reaction (I) with reactions (H), (IV), and (V). Although jkgOj^ was not 

measured in the work of reference 5, the change in HgO peak -AjllgOj relative to that ob- 
tained without added CO Ji^O^ was obtained. Inasmuch as 

Mo = k 4Mt 2 At + k 2[° H ]t[ H 2] t At <«> 

where |ohJ^ can be calculated by means of equation (5), tke P resence °f added 

CO is given by 



0 ) 


This value for jkgOj^ can tken * 3e usec * together with measured quantities to derive a sec- 
ond value kj for the rate constant of reaction (I). This derivation makes use of the fol- 


k l 

H 

t 

k 2 

h 9 
L 2 Jt 

; +k 4 

OH t 


( 10 ) 


And from equation (6), an expression for JohJ^. is obtained: 


[H = 


CO L 

cot 


1 


At 


(ID 
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From equations (10) and (11), the alternate value of the rate constant can be obtained: 



Equation (12) is an improvement over the equation used in reference 5 to calculate 
kp because all terms in equation (12) are known, or can be easily evaluated, whereas the 
previous equation contained a factor that was difficult to evaluate. 


RESULTS AND DISCUSSION 


Rate Constants for H 2 + OH Reaction 


Reaction-rate measurements of kg were made at or near each of the following tem- 
peratures: 304, 403, and 503 K. The experimental data and the rate constants derived 

from them are listed in table I. The value of k. needed for equations (3) and (5) was ob- 

12 ^ 

tained by setting it equal to 1. 55x10 cubic centimeters per mole per second at 300 K 
(ref. 2), with an assumed activation energy of 1 kilocalorie per mole (4.2 kJ/mole) 

(ref. 1). 

The rate constants listed in table I show a good deal of scatter. The reason for this 
is the great sensitivity of the results to the measured water concentration l^oj. Calcu- 
lations made using equations (2) and (5) and assuming 10 and 20 percent errors in HgC)J 
showed that the errors were magnified in the resulting kg. In general, 10 percent er- 
rors in [kgO] would produce 25 to 40 percent changes in kg; and 20 percent errors, 

40 to 60 percent changes. Furthermore, since the mass-spectrometer stability was 
about ±1/2 peak unit and HgO signals were usually only about 5 to 6 peak units, there was 
an inherent ±10 percent uncertainty in j^HgO] . hi table I the standard deviations show that 
the errors in HgO measurements must have averaged 10 to 15 percent, at most. The 
results are plotted in the Arrhenius fashion in figure 4, and the least-squares line through 
the data yields the following expression for the rate constant: 
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kg = (4. 7±1. 2)xl0 ^ exp 


~-(4200±300) 

RT 


cm'* mole * sec" 


(13) 


The error limits quoted are based on the standard deviations for the slope and intercept 
of our least-squares line. The activation energy of 4. 2±0. 3 kilocalories per mole 
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(17. 6±1. 3 kJ/mole) changes only slightly if other values are assumed for the activation 
energy of reaction (IV), as shown in the following table: 


Assumed activation 
energy, 
reaction (IV) 

Resulting activation 
energy, 
reaction (n) 

kcal/mole 

kJ/mole 

kcal/ mole 

kj/mole 

0 

0 

3. 5 

14.7 

. 5 

2.1 

3.7 

15. 5 

1.0 

4.2 

4.0 

16.7 

1. 5 

6.3 

4. 1 

17.2 

2.0 

8. 4 

4.3 

18. 0 


The effect of values on kg is also small. Analysis shows that a given percent- 
age change in values will produce only about one-half of that percentage change in kg. 

The dashed line in figure 4 represents the recommendation of the Leeds group 
(ref. 8). It is based on the room -temperature results of references 1, 2, and 9, plus 
eight other determinations made at higher temperatures in flame, shock-tube, and 
ignition- limit studies. All these high -temperature data lie above the extrapolation of the 
line derived from the present results, and the slope of the dashed line corresponds to an 
activation energy of 5. 15±0. 15 kilocalories per mole (21. 56±0. 63 kj/mole) (ref. 8). 

Although the weight of the evidence seems to favor the higher activation energy, there 
are two facts that support the lower value. First, the only other direct determination of 
kg, made in the temperature range from 813 to 843 K (ref. 10) and shown as triangles in 
figure 4, lie below the solid line. Second, in two shock-tube studies involving the com- 
bustion of hydrogen, observations were made that could not be explained by the "ac- 
cepted" high-temperature values of kg but instead called for lower values (refs. 11 and 
12). In reference 12, the authors commented that their experimental data supported a 
kg value of about 2. 8x10^ exp (-3895/RT) which the authors of reference 12 obtained 
from N. R. Greiner in a private communication. 

All but one of the eight high-temperature determinations used in reference 8 to gen- 
erate the dashed line of figure 4 are, in some sense, indirect. One was a shock-tube 
study in which rate constants for a large number of reactions were varied to obtain a fit 
to observed ignition delays, and the results were not very sensitive to the rate of reac- 
tion (II). In three other studies, the ratio of kg to another rate constant was obtained, 
so that the value of the second rate constant had to be supplied in order to derive kg. 

Three more used flames in which the disappearance of OH was followed, and it was as- 
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sumed that all Hg consumed was by reaction (n). The final determination also employed 
a flame but was more direct: free-radical concentrations were measured by electron 
spin resonance (ESR) spectrometry and water was determined by freezeout. However, 
the production of water by reaction (IV) was neglected, and the values derived for kg 
may therefore be too large. 

In appears that the discrepancy between the fitted activation energy (ref. 8) and the 
value measured directly in this work must be resolved by redetermination of kg at tem- 
peratures above 1000 K. 

Recalculation of Rate Constants for CO + OH Reaction 

Both methods of deducing the rate constant for reaction (I) from the experimental 
data require that a value for kg be supplied in order to calculate |ohJj. (eq. (7)) or 
^HgOj t (eqs. (8) and (9)). In equation (12), however, kg appears in a more direct man- 
ner. Thus, there could be a discrepancy between k^ and k' if the value assigned to 
kg is wrong. 

As stated in the INTRODUCTION, a discrepancy did occur in the data obtained above 
room temperature and was reported in reference 5. Recalculation of the CO + OH rate 
constant using equations (6) and (12) and the newly obtained kg values (eq. (13)) showed 
no discrepancy between k^ and kj values. This fact suggests immediately that equa- 
tion (13) may be more correct than kg = 2. 3x10 ^ exp(-5200/RT) (ref. 2); but a recalcu- 
lation of k!^ values using the latter kg and equation (12) from the present work also 
showed very good agreement in k^ and k^ values. This latter result must mean that 
the equation used to calculate k^ values in reference 5 was inadequate. 

It was noted in reference 5 that k^ decreased as the CO concentration in the stirred 
reactor increased. Recent sampling experiments involving CO and Hg showed that CO 
can produce this effect because CO reduced mass-spectrometer sensitivity to other mass 
peaks more than did Hg. Thus, the greater the CO concentration in the sample, the 
smaller the ion current for a given mass peak. To correct for this effect, the recalcu- 
lated kj and k'^ for each temperature were plotted against £coJ, as in reference 5; and 
a least-squares line was drawn through each set of data to obtain the rate constants at 
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J 


zero CO concentration. The results are as follow: 


Temperature, 

K 

Rate constant, 

cm^ mole"*' sec”* 

k l 

k i 

310 

5.6xl0 10 

5.9X10 10 

440 

10.1 

10. 7 

610 

13.6 

13.3 


These data form an excellent straight line in an Arrhenius diagram and are represented 
by the following expression: 


k 


1 


(3.5±l.l)xlO n exp 


'-(1120*80) '] cm 3 

RT J 


mole 


-1 


sec 


-1 


(14) 


Error limits quoted are standard deviations for the slope and intercept of the least - 
squares line. 

The activation energy of 1. 12 kilocalories per mole (4. 69 kJ/mole) does not differ 
greatly from the previously reported (ref. 5) value of 1. 2 kilocalories per mole (5. 02 
kj/mole), but is considered more reliable in view of the close agreement between k^ 
and kj that now exists. The only other reported, but still unpublished, direct meas- 
urement of the activation energy for reaction (I) yielded 0±1 kilocalorie per mole (0±4.2 
kj/mole) (ref. 13), in satisfactory agreement with the present value. A critical analysis 
of all data available on this reaction up to 1968 (ref. 14) led to a recommended activation 
energy of 1. 08 kilocalories per mole (4. 52 kj/mole), which differs negligibly from 1. 12. 
Thus, the temperature dependence of kj seems well established. 

However, reanalysis of the data has done nothing toward increasing the 300 K value 
of kj, which is only about 65 percent as big as several other direct determinations made 
at room temperature and summarized in reference 14. This difference is outside the ex- 
perimental error. A possible explanation suggested in reference 5 is that mixing is in- 
complete in the stirred reactor because of the speed of the reactions occurring. But 
upon reflection this explanation seems less likely. For a stirred reactor to work prop- 
erly, the mixing time must be considerably less than both the residence time and the 
characteristic reaction time of the most rapid reaction. Under the conditions of these 
experiments, even the random -walk diffusion time is considerably less than the reaction 
time and only one-tenth the residence time. In addition, the entering gas velocity is so 
high that an average gas particle makes more than 40 traverses before it exits from the 
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reactor, and efficient mixing must have been attained. Thus, it is hard to attribute our 
differences from the results of others to poor mixing, and the reason for the discrepancy 
remains unclear. 


caSSbi 


Note Added in Proof 


Recently, N. R. Greiner's work on the 


H 2 + OH 


-H 2 0 + H and CO + OH 


C0 2 + H reactions using a spectroscopic technique was published (ref. 15). For the tem- 
perature range of 300 to 500 K, rate constants kg and were determined to be 
(4. 06±1.26)xl0 12 exp ’-(4020±180 )/Rt] and (1.26±1. 18)xl0 U exp [-(230±140 )/Rt] cubic 
centimeters per mole per second, respectively. Greiner's kg value is in excellent 
agreement with that reported in the present investigation. The disagreement in k^ val- 
ue is not surprising since the kj value at 300 K in the present investigation was much 
lower than those reported by several other workers. 


SUMMARY OF RESULTS 

Rate constants for the reaction Hg + OH HgO + H were measured at tempera- 

tures of 304, 403, and 503 K. The room -temperature result was close to other values in 
the literature but the activation energy, measured directly for the first time, was 4. 2±0. 3 
kilocalories per mole (17. 6±1. 3 kj/mole), which was less than the Leeds group's esti- 
mate of 5. 15±0. 15 kilocalories per mole (21. 56±0. 63 kJ/mole) that was based on indi- 
rect evidence. Reanalysis of previous CO + OH »- COg + H data by using either the 

newly determined activation energy or the previous energy of 5.2 kilocalories per mole 
(21. 8 kj/mole) for kg brought the two sets of CO + OH rate constants into essentially 
perfect agreement and established the activation energy of this reaction as 1. 12 kilo- 
calories per mole (4. 69 kj/mole). The discrepancy in the previous high-temperature 

CO + OH ► COg + H data was due to the fact that the previous equation that was used 

to calculate kj was inadequate. 

The Arrhenius expressions obtained for the rates of these two reactions were as fol- 
lows: 
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H 2 + OH: k 2 = (4.7±1.2)X10 12 exp|J^M00)j cm 3 mole" 1 sec' 1 
CO + OH: = (3. 5±1. l)xlO U expBlH2M2ll cm 3 mole' 1 sec” 1 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, January 13, 1970, 

129-01. 
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